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We apply a multiple-time version of the reductive perturbation method to 



study long waves as governed by the Boussinesq model equation. By requir- 



ing the absence of secular producing terms in each order of the perturbative 



H ' scheme, we show that the solitary-wave of the Boussinesq equation can be 

written as a solitary-wave satisfying simultaneously all equations of the KdV 
hierarchy, each one in a different slow time variable. We also show that the 
conditions for eliminating the secularities are such that they make the pertur- 
bation theory compatible with the linear theory coming from the Boussinesq 
equation. 
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I. INTRODUCTION 



As is well known, the intermediate long-wave Boussinesq model equation 



Utt - Uxx + U. 



'XXXX 




XX 



0, 



(1) 



with u{x,t) a one-dimensional field and with the subscripts denoting partial differentia- 
tion, is completely integrable0 and has N-soliton solutions. In particular, its solitary-wave 
solution is of the formi 



where k is the wavenumber. 

To study the long waves of eq.(|ll), we are going to consider a perturbative scheme based 
on the reductive perturbation method of Taniuti,i modified by the introduction of an in- 
finite number of slow time-variables: r^, t^, ry, etc. Then, as a consequence of a natural 
compatibility condition, we have that any wave field satisfying the KdV equation in the 
time Ts must also satisfy all equations of the KdV hierarchy,^ each one in a different slow 
time variable. The main reason for introducing these time variables, as we are going to see, 
is that they allow for the construction of a perturbative scheme, valid for weak nonlinear 
dispersive systems, which is free of solitary-wave related secularities. 

In this paper, by making use of the perturbative scheme with multiple slow time-scales, 
we will show that the solitary-wave of the Boussinesq equation may be written, in the slow 
variables, as a solitary-wave solution to the whole set of equations of the KdV hierarchy, each 
one in a different slow time variable. This result follows both, from the general long-wave 
perturbation theory, and from the observation that the perturbative series truncates for a 
solitary-wave solution to the KdV hierarchy equations, rendering thus an exact solution for 
the Boussinesq equation. Furthermore, we will show that the conditions for the elimination 
of the secular producing terms make the perturbation theory compatible with the linear 
theory associated to the Boussinesq equation. 




(2) 



2 



II. THE MULTIPLE TIME EVOLUTION EQUATIONS 



The long-wave limit is achieved by putting 



k = en 



where e is a small parameter. Accordingly, we define a slow space 



e = e(x - t) 



as well as an infinity of slow time coordinates: 



Ts = eH ; Ts = eH ; ry = eh ; 



Consequently, we have that 



and 



d _ d 
dx 



d d n d ^ d J d 



dt dr-s (9r5 dry 

In addition, we make the expansion 



2 ^ 2 / 2 4 



(3) 



(4) 



(5) 



(6) 



(7) 



(8) 



and we suppose that U2n = U2n{C,, T3, T5, ...), n = 0, 1, 2, . . ., which corresponds to an extention 
in the sense of Sandri.! Substituting eqs.(@), (|^ and (|) into the Boussinesq equation (P, 
the resulting expression, up to terms of order e^, is: 
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{uof + 2e\oU2 + e^(2uoM4 + (^^2)^) + 



0. 



(9) 



We proceed now to an order-by-order analysis of this equation. At order e°, after an 
integration in ^, we get 
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as 



Quo 



duo d^uo 
9^ ~ 'd^ 



as = -- 



(10) 



which is the KdV equation. 
At order e^, eq.(0]) yields 



d_ 



(11) 



Using the KdV equation ([T0| ) to express dup/dr^, integrating once in ^ and assuming a 
vanishing integration constant, we obtain 

>dun 



du2 ^ 1 d'^U2 



dup Id^up 3 d^up 9dupd'^up 9, .o^'^p 



+ 



• (12) 



Equation ([l2D , as it stands, presents two problems. First, the evolution of Up in the time 
is not known a priori. The second problem is that the term (d^up/dC,^), as a source term, 
is a secular producing term when up is chosen to be a solitary-wave solution of the KdV 
equation. In the next sections we will be dealing with these two problems. 



III. THE KORTEWEG DE VRIES HIERARCHY 

As we have seen, the field up satisfies the KdV equation ( p!OD in the time t^. The 
evolution of the same field up in any of the higher-order times T2n+i can be obtained in the 
following way.i First, to have a well ordered perturbative scheme we impose that each one 
of the equations for Mor2„+i be e-independent when passing from the slow (Mo,^,T2n+i) to 
the laboratory coordinates {u,x,t). This step selects all possible terms to appear in Mot2„+i- 
For instance, the evolution of up in time T5 is restricted to be of the form 

Mors = a5M0(55) + Pb^pUp^^^ + (/^g + -f5)Uo^Uo^^ + S^ulup^ , (13) 

where as, (3^, 75 and ^5 are unknown constants. Then, by imposing the natural (in the 
multiple time formalism) compatibility condition 

(uOrs) = (Upr2n+l) > (14) 

^ ■' T2n + 1 ^ / T3 
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with Morg given by eq.(0), it is possible to determine any 'Uot2„+i, i-^- to determine all 
constants appearing in Mor2„+i • As it can be verifiedjl the resulting equations are those given 
by the KdV hierarchy. In particular, for mqts and Uorr we obtain respectively 



uom ~ - + 30{uofuo^ , (15) 



and 



M077 = a7[ - ^^0(75) + 14uoMo(50 + 42^0^^0(45) + lAO{uofuo^ 

+ TOMo^^tiQg^g - 28OM0M05M055 - TOiuo^f - 70{uoYuo^^^] , (16) 

where and aj are free parameters not determined by the algebraic system originated from 
eq.(|l^. These free parameters are related to different possible normalizations of the slow 
time variables. 



IV. HIGHER ORDER EVOLUTION EQUATIONS 



We return now to eq.([T2|) for U2. Substituting Uorg from eq. ([TSl) , we obtain 



+ 



1 

9 

-- + 20«5 



M0(5C) + 



2 + lOas 



(mo)^moc • 



(17) 



We see thus that the secular-producing term ^0(55) can be eliminated if we choose 05 = |. 
In this case, Eq. (|l7|) becomes 

(18) 



M2t3 + 3(Mo«2)c - ^U2^^i = -3(no)^Mo5 + MOMOCCC - Uq^Uo^^ 



From this point on, we are going to consider some specific solutions to our equations. 
First of all, we assume the solution of the KdV equation ([T0|) to be the solitary-wave solution 



uq = — 2/t^sech^ — AasK'^Ts + 9 



(19) 



where ^ is a phase. However, we have just seen that uq must satisfy also the higher order 
equation ([T5|) of the KdV hierarchy. Actually, as we are going to see, to obtain a perturbative 



scheme free of secularities at any higher order, we will assume that uq be a solitary-wave 
solution to all equations of the KdV hierarchy, each one in a different slow-time variable. 
Such a solution is given by 

Uq = — 2K^sech^ — Aa^K'^T^ + IGa^tC'r^ — QAa-jti^r-j + ■ • • . (20) 

Using this solitary-wave solution, we see that the right-hand side of eq. ([18|) vanishes, leading 
to 

U2t3 + 3(moM2)5 - ^M2555 = , (21) 

which is a homogeneous linearized KdV equation. We will assume for it the trivial solution 

M2 = . (22) 
At order e^, and already assuming that M2 = 0, eq.(||) gives 

M4r3€ + 3(Mo«4)55 - ^M4(4g) = -Morr? + ^Orgrs • (23) 

Using equations (|TUp and (|T5|) to express Morg and Worg respectively, and integrating once in 
^, we obtain 

, , 1 1 45, , 2 

z io o 

35 5 75 45 15 

- y"05€"0555 - 2"o5"o{40 + -^uqUq^uq^^ - —{uofuo^ + — (mo?)^ • (24) 

The source term proportional tto{70 is the only resonant, that is, secular producing term to 
the solution u^. Then, in the very same way we did before, we first use the KdV hierarchy 
equation ( pj]) to express uorr After we do that, we can then choose the free parameter a-j in 
such a way to eliminate the resonant term from the right-hand side of eq. (^^ . This choice 
corresponds to ^7 = — ^, which brings eq.(^) to the form 

M4t3 + 3(MoM4)c - = ■^[^05MO{40 - M0«0(5O 

Z o 

+ IOuqUo^uq^^ - 5{uo^Y + 10(uo)^no^^^ - 20{uofuo^] . (25) 



Substituting again the solitary- wave solution (^) for uq, we can easily see that the nonho- 
mogeneous term of eq.(^) vanishes, leading to 

M4r3 + 3(MoU4)g - ^M45?C = ^ • (26) 

And again, we take the trivial solution 

M4 = . (27) 

It is easy to see that this is a general result that will repeat at any higher order: for n > 1, 
the evolution of U2n in the time rs, after using the KdV hierarchy equation to express "Uora^+i 
and substituting the solitary- wave solution ( pOf ) for uq, is given by a homogeneous linearized 
KdV equation. Consequently, the solution U2n = 0, for n > 1, can be assumed for any higher 
order. 



V. BACK TO THE LABORATORY COORDINATES 



Let us now take the solitary-wave solution to all equations of the KdV hierarchy. 



Uq = -2K^sech2 f< + 2k^T3 + 2k^t^ + Ak't-j + 



(28) 



where we have already substituted the (not anymore) free parameters a2n+i, and rewrite it 
in the laboratory coordinates. First, recall that we have expanded u according to eq.(§). 
Thereafter, we have found a particular solution in which U2n = 0, for n > 1. Consequently, 
expansion (H) truncates, leading to an exact solution of the form 



U = t Uq 



(29) 



with Mo given by eg . (pS]) . Moreover, the slow variables (/t, ^, T2n+i) are related to the labora- 
tory ones, {k,x,t), respectively by eqs.(|^), @ and (^. Then, in the laboratory coordinates, 
the exact solution f^Ol) is written as 



u = — 2/c^sech^ k 



X 



- (l - 2k^ - 2k^ - Ak^ 
7 



(30) 



Now, the series appearing inside the parenthesis can be summed: 

l_2k''-2k^-Ak^---- = (l- A^y^' . (31) 

Consequently, we get 

u = -2fc2sech2 [k (a; - Vl - Ak^ t)] , (32) 
which is the well known solitary-wave solution of the Boussinesq equation (|l|). 

VI. RELATION TO THE DISPERSION RELATION EXPANSION 

Let us take now the linear Boussinesq dispersion relation: 

u{k) = k {l + k^y^\ (33) 
Its long-wave {k = en) expansion is given by 

uj{k) = eK + a^e'tC' + a^e^ti^ + aye^fi;^ H , (34) 

where the coefficients a2n+i, except for which arose naturally in the KdV equation (|10|), 



coincide exactly with those necessary to eliminate the secular producing terms in each order 
of the perturbative scheme. With this expansion, the solution of the associated linear 
Boussinesq equation can be written as 



u = exp I 



Ke{x -t) + a-sK^e^t + a^n^e^t + aTK^e^t H . (35) 



Therefore, if we define from the very begining, as given by this solution, the properly nor- 
malized slow time coordinates 

T3 = a^eH ; n = a^e^t ; T7 = a-je't ; . . . , (36) 

the resulting perturbative theory will be automatically free of secularities.0 Furthermore, 
the linear limit of the perturbation theory will be compatible with the linear theory coming 
directly from the Boussinesq equation (|l|). 
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VII. FINAL COMMENTS 



By applying a multiple time version of the reductive perturbation method to the Boussi- 
nesq model equation, and by eliminating the solitary-wave related secular producing terms 
through the use of the KdV hierarchy equations, we have shown that the solitary-wave of 
the Boussinesq equation is given by a solitary-wave satisfying, in the slow variables, all 
equations of the KdV hierarchy. Accordingly, while the KdV solitary-wave depends only 
on one slow time variable, namely t^, the solitary-wave of the Boussinesq equation can be 
thought of as depending on the infinite slow time variables. 

The same results hold! for the shallow water wave (SWW) equation. In other words, the 
solitary-wave of the SWW equation can also be written as a solitary-wave satisfying simulta- 
neously, in the slow variables, all equations of the KdV hierarchy. It is important to remark 
that in both cases the resulting secularity-free perturbation theory will be automatically 
compatible with the corresponding linear theory. 

A crucial point of the multiple time perturbative scheme is the return to the laboratory 
coordinates, which implies in a renormalization of the solitary-wave velocity.^ In the case 
of the Boussinesq, as well as of the SWW equation, this renormalization is such that the 
KdV hierarchy solitary-wave is led to the corresponding Boussinesq or SWW solitary-waves. 
However, when the original nonlinear dispersive system does not present an exact solitary- 
wave solution, the series will not truncate. In this case, a secularity-free expansion can still 
be obtained and the process of returning to the laboratory coordinates can be made order- 
by-order at any higher order, implying in a sucessive renormalization in the velocity of the 
solitary-wave, which in this case is represented by the leading order term of the perturbative 
series. i 
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